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1.0  INTRODUCTION 


Nuclear  magnetic  resonance  (NMR)  gyro  developments  over  the  last  20 
years  have  been  aimed  at  providing  a  low  cost,  high  reliability 
alternative  to  conventional  mechanical  gyroscopes.  By  removing  the  need 
for  precision  moving  parts  such  devices  would  be  inherently  insensitive  to 
mechanical  shock  and  vibration;  and  thus  well  suited  to  strapdown 
systems.  Two  major  NMR  gyro  developments  in  the  United  States  produced 
engineering  models  which  demonstrated  rate  bias  stabilities  of  better  than 
1  deg/hr  '  .  The  goals  for  these  devices  ranged  from  0.1  deg/hr  to  0.01 
deg/hr.  However,  the  success  of  the  laser  gyro,  which  was  aimed  at  the 
same  market,  caused  abandonment  of  both  of  these  NMR  gyro  programs. 

Por  high  accuracy  applications,  it  has  long  been  recognized  that 
conventional  inertial  navigation  system  technology  has  advanced  to  the 
point  where  uncertainties  in  the  knowledge  of  the  earth’s  shape  and 
gravity  field  represent  significant  sources  of  navigational  error3. 
Thus,  efforts  are  underway  to  develop  gravity  gradiometers  for 
incorporation  in  high  accuracy  inertial  measurement  units.  A  room 
temperature  gravity  gradiometer  has  completed  its  initial  sea  trials  and 
demonstrated  an  improved  on-line  measurement  of  the  gravity  gradient. 
This  improved  measure  of  gravitational  acceleration  will  reduce  the  error 
for  the  deduced  inertial  acceleration;  and  the  next  generation  inertial 
measurement  unit,  incorporating  such  a  device,  may  well  provide  improved 
resolution  using  room  temperature  technology.  However,  potential  advances 
possible  with  cryogenic  instruments  may  well  exceed  any  room  temperature 
technology  in  the  years  ahead. 


To  realize  the  advantage  possible  with  any  one  cryogenic  instrument 

requires  development  of  an  entire  family  of  cryogenic  inertial 

instruments,  since  all  must  be  integrated  onto  a  single  isothermal 

platform.  Stanford  University  is  developing  the  technology  base  to 

support  both  cryogenic  gravity  gradiometers  and  cryogenic  accelerometers. 

The  3He  nuclear  gyro,  the  subject  of  this  report,  is  the  candidate  gyro 

for  such  an  all  cryogenic  inertial  measurement  unit.  Stanford  University 

3 

has  recently  designed  and  constructed  an  all  fused  quartz  He  gyro 
housing  consistent  with  a  single  axis  angular  stability  approaching 
2x10  5  deg/hr  (Figure  1-1). 


2.0  BASIC  PRINCIPLES 


3 

The  He  nucleus  possesses  both  intrinsic  spin  angular  momentum  and  a 
magnetic  dipole  moment  which  is  directed  antiparallel  to  the  spin  axis. 
However,  a  sample  of  3He  will  not  generally  possess  any  net  angular 

momentum  or  magnetization  due  to  the  random  orientation  of  the  individual 

4 

spins.  The  process  of  optical  pumping  is  employed  to  orient  the 
individual  spins  along  a  preferred  direction  and  thus  achieve  a  net  sample 
angular  momentum. 

When  such  a  polarized  sample  is  placed  in  a  uniform  magnetic  field,  B, 

two  processes  ensue.  The  first  is  a  relaxation  of  the  sample  back  to  its 

unpolarized  equilibrium  condition  (assuming  the  optical  pumping  process  is 

terminated).  The  equilibrium  value,  V  ,  is  determined  using  Boltzmann 

o 

statistics  to  be 


N |y | tank 


(2-1) 


3 


I 


Where  N  is  the  number  of  3He  atoms,  y  is  the  3He  magnetic  dipode 

moment,  k  is  Boltzmann’s  constant  and  T  is  the  temperature.  At 
B 

equilibrium  the  net  polarization  lies  along  the  applied  field  direction. 
Given  that  the  equilibrium  condition  is  approached  at  a  rate  proportional 
to  the  displacement  from  equilibrium,  the  relaxation  process  can  be 
described  by  the  differential  equation 
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where  it  is  assumed  that  the  applied  field  is  in  the  z  direction.  The 
characteristic  time  Tj  is  referred  to  as  the  longitudinal  relaxation 
time  and  T 2  is  the  transverse  relaxation  time.  T  and  T 2  may  be 
nearly  equal;  but,  in  general  T^  is  greater  because  of  magnetic  field 
gradient  effects. 

The  second  process  is  due  to  the  applied  field  interacting  with  the 
sample  magnetization  to  produce  a  torque,  M  x  B.  This  torque  equates  to 
the  time  rate  of  change  of  the  sample  angular  momentum  in  an  inertial 
frame; 


d 

dt 


H 


M  x  B 


(2-3) 


But,  H  and  M  are  antiparallel  and  in  fact  are  related  by  the  gyromagnetic 
ratio,  Y,  which  is  a  constant  for  the  species. 


M  =  YH 


(2-4) 


Thus  the  equation  of  motion  can  be  written  solely  in  terms  of  M: 


It  x  yt 


(2-5) 


This  equation  holds  for  an  inertial  frame:  but  T?  is  tied  to  the  gyro. 
Then  if  the  rotation  rate  of  the  gyro,  with  respect  to  the  inertial  frame, 
is  co  ,  the  equation  of  motion  in  the  gyro  frame  is: 


Mx  yB  - 


<o  x 


(2-6) 


=  M  x  (  y ^  +  <o  ) 


(2-7) 


Then  combining  equations  (2)  and  (7),  with  "b  still  assumed  to  lie  in  the  z 
direction,  yields 
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3  4 

For  a  liquid  sample  of  He  in  He,  a  T  of  greater  than  140  hours 

has  been  obtained5.  In  principle  T  can  be  made  on  the  order  of 

weeks.  Thus  for  many  calculations  the  relaxation  effects  are  neglected 
leaving 


(2-9) 


This  equation  is  the  basis  of  the  nuclear  gyro.  The  solution  indicates 

2  2 

that  M  precesses  about  an  axis  -(YB-Kd)  at  a  rate  of  [(YB+w^)  +^x+ 

2  1/2 

(o  y]  rad/sec.  Since  Y  is  a  constant  for  a  given  nuclear 

species,  the  Larmor  precession  frequency,  y  B,  can  be  predicted.  Any 
deviation  from  the  Larmor  rate  is  assumed  to  be  due  to  a  non-zero  ^  . 

Equation  (9)  shows  that  the  greatest  sensitivity  of  the  gyro  occurs  for 

A 

rotations  in  the  2-direction;  but  it  is  also  clear  that  there  is  a 
mechanism  to  sense  cross  axis  rates  as  well. 

In  order  to  utilize  these  principles  in  a  practical  device  requires 
several  elements.  The  first  is  a  nuclear  species,  with  intrinsic  spin 
angular  momentum.  Ideany  the  species  should  be  a  spin  1/2  species; 
otherwise  it  will  possess  an  electric  quadrupole  moment.  An  electric 
quadrupole  moment  is  due  to  the  ellipticity  of  the  distribution  of  charge 
in  the  nucleus.  The  important  point  is  that  such  a  moment  will  interact 


with  an  electric  field  gradient  to  produce  a  torque  on  the  nucleus, 
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leading  to  a  shift  of  the  tarmor  frequency.  A  spin  1/2  species,  however, 

has  only  a  magnetic  dipole  moment}  and  thus  the  only  interaction  of 

concern  will  be  with  a  magnetic  field.  3He  is  such  a  species.  It  also 

has  the  added  advantage  of  being  useable  as  either  a  liquid  or  a  gas  at 

cryogenic  temperatures,  and  demonstrates  very  long  relaxation  times. 

The  second  element  required  is  a  means  of  polarizing  the  nuclear 

3 

species.  Optical  pumping  has  been  successfully  employed  on  He  at  room 
4  3 

temperature  .  A  gas  He  sample,  polarized  at  room  temperature,  can  be 

4 

condensed  in  solution  with  He  without  significant  loss  of 
.  5 

polarization  . 

The  third  element  is  an  applied  magnetic  field.  One  obvious 

difficulty  in  making  a  practical  3He  nuclear  gyro  is  in  obtaining  a 

uniform,  static  magnetic  field.  The  magnitude  of  the  earth's  magnetic 

-5 

field  is  on  the  order  of  0.5G  (5x10  T)  at  the  surface,  and  varies  in 

magnitude  and  direction  from  one  location  to  another.  So  one  task  is  to 

screen  the  3He  sample  from  ambient  magnetic  fields.  A  method  has  been 

—8  —12 

developed  to  create  low  field,  10  G  (10  T),  regions  within  a 
superconducting  lead  foil  shield6.  This  type  of  shield  has  the  property 
of  perfect  diamagnetism.  That  is,  beyond  a  penetration  depth,  which  is 
small  compared  to  the  foil  thickness,  the  magnetic  flux  density  in  the 
foil  is  zero.  Thus  changing  magnetic  field  environments  external  to  the 
foil  have  absolutely  no  effect  within  the  shielded  region. 

In  addition  to  shielding  from  outside  magnetic  fields,  a  high  degree 
of  uniformity  is  required  of  the  applied  magnetic  field.  If  sixth  order 
superconducting  Helmholtz  coils  are  used  to  generate  the  applied  field, 
the  field  gradiants  will  be  largely  due  to  inhomogeneities  in  the  field 


7 


trapped  within  the  lead  foil  shield.  Based  on  current  capabilities  for 

establishing  ultra-low  field  regions,  the  major  field  inhomogeneity  is 

—8  —1 

estimated  to  be  a  linear  gradient  of  magnitude  less  than  10  G-cm 

(10  10T-m  Figure  2-1  shows  a  three  axis  magnetic  field  profile 

of  an  ultra-low  field  trapped  within  a  20  cm  diameter  lead  foil 
6 

shield  .  Furthermore,  the  field  generating  coils  can  be  made 
superconducting  so  that  the  applied  field  will  be  generated  a  persistent 
supercurrent,  which  is  inherently  very  stable. 

The  remaining  element  of  the  ^He  nuclear  gyro  is  the  readout 
mechanism.  A  sensitive  magnetometer  is  required  to  monitor  the  precessing 
magnetization  vector.  Since  there  is  already  a  requirement  to  maintain 
the  shield  and  the  field  coils  at  cryogenic  temperatures,  a  natural 
candidate  for  the  magnetometer  is  the  SQUID  (Superconducting  Quantum 
Interference  Device)  magnetometer.  Fortunately  this  is  also  the  most 
sensitive  magnetometer  currently  available.  Theory  and  use  of  the  SQUID 
are  described  in  a  number  of  references  '  '  '  s  SQUID  use  for  a  He 
nuclear  gyro  has  been  described  by  Taber5. 


S 


Figure  2-1.  Three  Mutually  Perpendicular  Magnetic  Field  Component 
Profiles  Along  the  Axis  of  a  20  cm  Diameter 
Super-conducting  Shied 


3.0  CROSS  COUPLING  EFFECTS 

It  is  clear  from  Equation  (2-9)  that  the  nuclear  gyro  is  not  a  true 
single  degree  of  freedom  gyro.  Indeed,  while  the  axis  along  the  applied 
field  has  the  greatest  sensitivity,  the  device  is  affected  by  cross  axis 


rates 

as 

well11.  So, 

while 

we 

are  trying  to  measure  co^,  the 

instrument 

is  sensitive 

to  0) 

X 

and 

to  also.  Take,  for  example,  the 

case 

where 

0)  =  0)  = 

x  z 

0  and 

0) 

y 

0  then  solving  Equation  (2-9) 

yields  the  following 


M^(  t )  cos  cot 


^  sin  (ijt 

£0 


y  sin  cot 
0) 


vo) 


M  (t)  «  -  -^sin  cot  cos  cot  +\co  )  (l-cos  cot)  — -■£  (1-cos  cot)  M  (0)  (3-1) 

y  co  &)2  y 

<u  YBco  ,  rX  2 

M  (t)  — ¥•  sin  cot  - -2  (1-cos  cot)  cos  cot  (1-cosco  t)  M  (0) 

z  co  2  'co  /  — 1 — z  — 1 

— J  l—  co 


where 


co  =  Y  (yB)2  +  co2 


One  effect  of  a  cross  axis  input  is  immediately  apparent.  Even  with  no 

component  of  rotation  about  the  applied  field  (input)  axis,  the  precession 

frequency  is  y<YB):  "+  co2  rad/sec.  This  must  be  accounted  for  somehow  to 

2 

avoid  a  gyro  drift  error  on  the  order  of  1/2  (o>  ^/yB)  rad/sec. 

Figure  3-1  shows  the  effect  of  an  uncorrected  cross  axis  input.  As 

2 

shown  above  the  drift  grows  as  cu^  and  inversely  with  the  Larmor 
frequency,  y  B.  This  suggests  that  cross  coupling  effects  can  be  made 


10 


arbitrarily  small  simply  by  increasing  B.  However,  for  a  cross  axis  input 

rate  of  1.0  rad/sec,  the  Larmor  frequency  must  be  107  rad/sec  to  reduce 

—8  3 

the  uncorrected  drift  to  5  x  10  rad/sec  (0.01  deg/h).  For  He  this 

-2 

requires  that  B  be  about  490  G  (4.9  x  10  T).  The  gradients  associated 

with  such  large  fields  make  them  impractical.  The  better  approach  is  to 
measure  the  cross  axis  inputs  and  then  compensate  for  their  effects. 

So  consider  the  general  case  where  there  is  both  an  input  axis  and  a 
cross  axis  component  of  the  input  rate.  Take  the  two  respective  components 


of  &  to  be  to  and  to  , 
z  y 


There  is  no  loss  of  generality  here  since  the 


direction  of  the  Y-axis,  in  the  plane  normal  to  the  Z-axis,  is  arbitrary. 


Proceeding  as  before,  we  get 


d_ 

dt 


M 


The  solution  is  thus  the  same  as  (3.1)  with  the  term  yB  replaced  by  (y  B 

+  ^M.  Thus  the  precession  frequency  seen  in  the  gyro  frame  is  now  y[(B 
2  2  1/2 

+  to  /  y  )  +  (  to  /  y )  ]  '  .  it  can  also  be  shown  that  the  normal 

z  y 

to  the  tip  path  plane  (precession  plane)  of  M  is  parallel  to  y  (B  +  to^/y ) 

■A 

Z  +  y  (  to^/y  )  Y.  That  is  the  behavior  which  is  identical  to  the  case 

where  w  =  0  and  the  applied  field  is 


B 


to  .  to 

(B  ♦  y-  >  Z  ♦  (-^) 


A 

Y. 


11 


7B  (rad/sec) 


Figure  3-1.  Drift  Rate  vs.  barmor  Frequency  for  a  Number  of  Cross 
Axis  Input  Rates,  oy. 


where  Z  and  Y  are  unit  vectors  along  the  Z  and  Y  axes  respectively.  This  is 

simply  a  statement  of  the  equivalence  between  rotations  and  magnetic 
12 

fields  .  Thus  the  effect  of  a  cross  axis  rotation,  &)y,  is  the  same  as 
that  of  applying  a  magnetic  field,  o^/Y  »  in  the  same  direction  as  the 

«*■  -A  A  A  A 

cross  axis  component  of  <»  .  Thus,  if  u  =  w  x  +  W  Y  +  w  z,  the 

x  y  z 

sensed  precession  frequency,  ft,  will  be  given  by 


222  1/2 
ft  =  r  (Y  8  +  +  oj  +w]  . 

z  x  y 


(3-3) 


SIGNAL  LOSS  DUE  TO  CROSS  AXIS  INPUTS 


Another  effect  of  a  cross  axis  rotation  is  to  tilt  the  precession  plane 


such  that  it  is  normal  to  the  effective  field,  B 


0)  J.  U)  A  0)  . 

B  =*-£  X  +  -2  y  +  (B  +  -£)  Z. 
eff  Y  Y  Y 


(3-4) 


The  instantaneous  normal  to  the  precession  plane  will  always  be  in  the 

direction  of  B  however,  the  manner  in  which  the  normal  tracks  B  __ 

eft  eft 

depends  on  how  fast  Bef£  changes  direction  and  on  the  inertial  conditions 

^  —a 

of  M  when  the  change  transpires.  If  M  initially  lies  in  the  plane  of 
precession,  this  condition  can  be  maintained  provided  the  rate  of  change  of 


direction  of  B  ^  is  small  compared  to  |YBeff  I  •  Violation  of  this 
condition  will  result  in  M  tracing  out  a  cone  around  B  .  As  an  example, 
consider  the  extreme  case  where  initially  B  ^  *  B  Z ,  then  as  M  points 

-A  A 

along  the  Y-axis,  there  is  an  instantaneous  change  in  w  to  be  Y.  This 

a  _L 

tilts  the  effective  field  toward  the  Y  direction.  Now  M  maintains  a  fixed 


n 


angle  with  respect  to  B  __  and  thus  moves  up  out  of  the  horizontal  plane 

er  r 

(see  Figure  3-2).  If  at  some  time  later,  when  M  is  at  its  maximum  height, (u 
goes  back  to  zero,  the  precession  plane  will  again  be  a  horizontal  plane; 

however,  M  no  longer  lies  in  this  plane.  In  fact,  if  the  instantaneous 

_ _ 

change  in  co  rotated  B  by  the  angle  0  ,  then  in  the  final  condition  M 

er  t 

will  make  an  angle  20  with  the  horizontal  plane.  Recall  that  the 
magnetometer  senses  the  horizontal  component  of  M.  Thus  for  the  example 
above  the  signal  has  been  reduced  to  have  a  peak  amplitude  of  | "fi |  cos  20  . 

A 

z 


Figure  3-  2.  Signal  Reduction  Due  to  Cross  Axis  Rate. 

Initially  co  =  0,  then  at  tj,  (o  =  (Uy  Y. 

This  tilts  the  precession  plane.  At  time 
1 2 » again  becomes  zero  and  the  precession 
plane  becomes  horizontal.  Note  now  that  M  has 
a  constant  Z  component;  and  the  X  component  has 
been  reduced. 
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EPFECT  OF  OSCILLATORY  CROSS-AXIS  RATE  INPUTS 


A  problem  that  has  arisen  involves  oscillatory  rates.  In  particular,  if 
the  nuclear  gyro  experiences  an  oscillatory  rate*3  about  a  cross  axis  and 
at  a  frequency  corresponding  to  the  Larmor  frequency,  a  periodic  loss  of 
magnetometer  output  results.  This  was  demonstrated  with  a  digital 
simulation  of  the  3He  gyro  dynamics.  The  particular  case  demonstrated  was 
an  X-axis  gyro  in  response  to  an  oscillating  rate  about  the  Y  axis.  The 
cross-axis  rate  input  (Figure  3-3)  is  described  by 


0) 

y 


0.05  yB  sgnfcos  YBt] 


(3-5) 


The  resulting  path  of  the  magnetization  vector,  as  seen  in  the  gyro  frame, 
is  shown  to  wind  its  way  up  the  gyro  X  axis  (Figure  3-4).  The  magnitude  of 
M  is  unchanged  since  relaxation  effects  are  ignored,  but  if  the  oscillatory 
input  persists,  M  winds  back  down  into  the  X-0  plane  and  then  continues  its 
way  down  the  -X  axis  and  back  and  forth.  Since  the  magnetometer  senses  the 
component  of  M  orthogonal  to  B  the  effect  is  to  alternately  diminish  and 
restore  the  magnetometer  signal.  Overall  then,  the  signal-to-noise  ratio  of 
the  instrument  is  decreased.  In  this  demonstration  time  has  been  scaled  in 


terms  of  yB. 


For  example,  if  YB  is  chosen  to  be  2  tt  rad/s. 


then  to  has  a 

y 


peak  rate  of  0.31  rad/s  and  the  total  time  shown  (Figure  3-4)  is  about  7.5 


s.  In  other  words,  if  the  peak  input  rate  is  5%  of  the  Larmor  rate,  it  will 


take  about  7.5  Larmor  periods  for  M  to  wind  its  way  out  of  the  X*0  plane. 
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The  situation  is  most  easily  explained  in  a  rotating  frame.  Por 
example,  approximate  w  as 

u  =»  0.05  YB  cos  YBt  (3-6) 

y 

Then  represent  (a  as  two  counter-rotating  vectors  in  the  X=0  plane,  each 
of  magnitude  0.25  YB,  starting  in  the  Y-axis  direction  and  rotating  at  YB 
rad/s.  Then,  in  the  rotating  frame,  the  vector  traveling  with  the  rotating 
frame  is  stationary,  while  the  other  appears  to  rotate  at  twice  the  Larmor 
rate.  The  effects  of  the  latter  average  to  zero,  but  the  former  causes  M  to 
precess  about  the  rotating  frame  Y  axis.  In  the  gyro  frame  both  the  Y  axis 
motion  and  the  Larmor  precession  about  the  X  axis  cause  the  spiral 
trajectory  shown.  This  situation  is  a  worst-case  situation.  For 
oscillatory  rates  much  higher  or  lower  than  YB»  the  effect  averages  to  zero. 


Figure  3-3.  Simulated  Oscillatory  Y-axis  Rate 


Figure  3-4.  Effect  on  X-axis  Gyro  of  Oscillatory  Rate  about  Y-axis. 

RESOLUTION  WITH  THREE  ORTHOGONAL  GYROS 

In  general,  the  uncorrected  drift,  [ft  -  (  Y  B  +  oi  )]  is  not  tolerable. 

z 

So  consider  using  three  identical  gyros,  rigidly  aligned,  such  that  their 
input  axes  are  orthogonal  (Figure  3-5).  Then  each  gyro  will  undergo  the 
same  rotation  but,  because  the  input  axes  are  orthogonal,  each  will  be 
affected  differently. 

_  A  A  A 

For  the  case  where  oj  =  X  +  Y  +  2,  the  precession 

frequencies  sensed  by  each  gyro  are 
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fc 

k 

k 


I 


ft  2  =  (Y  B  +  W  )2  +  O)2  +  CO2 
x  x  y  z 


( 3-7a) 


ft  2  =  w2  +  (Y  B  +  w  )2  +  W2 
y  X  y  z 


( 3-7b) 


ft  2  =*  (o2  +  a)2  +  (  yb  +  co  ) 2 

z  x  y  2 


(3-7 c) 


It  is  convenient  to  normalize  by  dividing  both  sides  of  the  above  by 
(  Y3)2,  giving: 


ft 


ft 


ft 


2 

x 


2 

y 


2 

z 


u  +v 

a)2  +  (1 
— x 


2 


2  2 
+  (0  +.  (1) 

— y  -- z 


2  2 
0)  )  +6J 

-y  -  z 


+  (1  +  OJ  ) 

— z 


2 


( 3-8a) 

( 3-8b) 

( 3-8c ) 


where  the  underbars  indicate  normalized  quantities. 


]  8 


<x 


Now  the  question  is,  do  Equations  (3-8)  uniquely  determine  "5  ,  and,  in 
general,  the  answer  is  no.  Knowledge  of  the  time  history  of  to  will  not 
resolve  the  ambiguity,  as  the  following  example  will  show.  Suppose 

-  A  A  ^ 

initially  (o  =  -1/3  X  -  1/3  Y  -  1/3  Z.  This  particular  case  is  resolved 
unambiguously  from  (3-8).  Furthermore,  assume  that  this  condition  persists 
for  some  time.  Then  suppose  the  input  changes  such  that  =  -0.3333  X  - 

A  A 

0.3333  Y  -  0.3333  Z.  With  only  the  measurements  f)  ,  ft  ,  Q  and 

—  x  —  y  —  z 

knowledge  of  y  B,  two  solutions  for  are  possible.  One  solution  is  indeed 

_ k  A  A 

the  true  solution;  the  other,  however,  is  u  =  -0.  33337  X  -  0.33337  Y  - 

A 

0.33337  Z.  The  two  solutions  are  nearly  equal,  though  in  opposite 
directions  from  the  previous  known  solution.  The  smaller  the  step  away  from 

-A  A 

-  1.3  X  -  1/3  Y  -  1/3  Z  the  closer  are  the  two  solutions.  Thus  it  is  easy 
to  conceive  of  the  possibility  of  locking  onto  the  false  solution  from  this 
starting  point. 

Let  us  investigate  the  source  of  the  ambiguity  by  considering  the  space 

of  possible  inputs,  Sb  1,  0)  y »  —  z'*’  in  which  we  graph  Equation 

(3-8).  Figure  3-6  shows  the  case  where  go  =*  0.  Here  the  false  solution  is 

(-2/3,  -2/3,  -2/3),  and  =  2.  y  =  =  1.0.  We  see  that  the 

information  from  any  one  gyro  defines  a  sphere  on  which  the  possible  values 

of  us  lie.  The  solution(s)  then  will  be  the  point(s)  common  to  all  three 

spheres.  Figure  3-6  is  an  isometric  drawing  and  so  the  two  solutions  plot 

as  the  same  point;  nevertheless  one  can  easily  visualize  the  false  solution 

for  this  care.  Note  that  for  three  identical  gyros  that  the  centers  for  the 

three  spheres  are  located  at  (-1,  0,  0),  (0,  -1,  0),  and  (0,  0,  -1)  and  the 

corresponding  radii  are  ft  ,  12  ,  and  ft  . 

—  x  —  y  — z 
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What  then  is  the  nature  of  the  ambiguity?  Consider  the  case  where  (o 
lies  outside  the  plane  of  centers.  Choosing  any  two  of  the  spheres  we  see 
that  their  intersection  must  be  a  circle.  Furthermore,  the  plane  of  centers 

cuts  this  circle  in  half;  i.e.,  the  center  of  this  circle  lies  in  the  plane 

* 

of  centers.  Now  since  the  third  sphere  must  also  pass  through  co  ,  this 
sphere  will  intersect  the  circle  at  two  points.  One  is  indeed  at  Sf  and  the 
second  point  is  the  reflection  of  cu  across  the  plane  of  centers. 

We  can  now  propose  a  simple  criterion  for  resolving  ambiguities.  Since 
one  would  want  to  stay  away  from  the  resonance  condition  ( w  =  y  b)  in 
general,  the  Larmor  frequency  could  be  chosen  such  that 

0)  +  a)  +  to  >  -yB  (3-10) 

x  y  z  — 

for  all  allowable  of.  Then  the  proper  choice  of  solutions  will  always  be 
on  the  same  side  of  the  plane  of  centers.  Another  possible  criterion, 
though  more  restrictive,  would  be  to  choose 

Iy  B  j  >  ySj^J  max  ;  (3-11) 

then  the  proper  solution  will  be  the  one  with  the  magnitude  less  than  1/V  3 
(  |  y  B  j  ). 

USE  OF  ADDITIONAL  MAGNETIC  FIELDS  TO  NULL 
CROSS  AXIS  INPUTS 

It  is  clear  from  Figure  3-4  that  some  strategy  must  be  employed  to  keep 
the  magnetization  in  the  precession  plane  despite  the  presence  of  cross  axis 
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rotations.  We  could  consider  restricting  the  operating  regime  to  be  always 
far  away  from  resonance  conditions;  but  there  may  be  yet  a  more  satisfactory 
approach.  Consider  again  the  Bloch  equation  for  an  arbitrary  3 . 


A  A 

If  we  add  coils  for  generating  magnetic  fields  in  the  X  and  Y  axes  we  get 


since  the  nuclear  gyro  is  sensitive  to  cross  axis  inputs,  it  is 
theoretically  possible  to  determine  all  three  components  of  a?  with  a 
single  device.  To  do  this  three  magnetometers  are  required,  one  to  measure 
each  components  of  M.  The  procedure  for  determining  3  from  measurements 
of  the  components  of  M  is  as  follows: 

1)  Take  three  consecutive  readings  of  "m,  m  ,  m!  ,  and  if" . 

a  b  c 

2)  Form  the  vectors  "m,  ■  M,  -  M  and  $  *  M  -if. 

Aba  B  c  b 
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3)  Assuming  6J  is  constant  during  the  three  readings,  then  the  vectors 

M  and  "m  lie  in  the  precession  plane  and  x  M0/)  M^  j  [Mg]  is  a 

A  A. 

unit  vector  in  the  direction  -(  YBy  +  x  B  +  <*>  )  *  *  (Yb 

X  A  2  J 

A  A 

+  0) )  Z  which  is  denoted  as  n. 

2 

_ fc 

4)  Next  choose  any  of  the  three  measurements,  say  Mb,  and  determine 

m!  x  n  =  I  M.  I  sin  <J>  or 
b  1  b 


(3-14) 


5)  The  radius  of  the  precession  plane  is  given  by  |  Hb  |  sin  <(>  = 

fib  x  "  I  * 

6)  Now  [ co |  can  be  determined  (see  Figure  3-8),  from  the  radius  of  the 
precession  plane,  the  time  between  measurements,  and  the  measurements  Mg 


Figure  3-8.  Determination  of  0)  from  three  measurements  of  M. 

The  veiw  is  directly  into  the  precession^plane. 

Note  At  is  the  time  between  measurements  Mc  and  Ma. 

Assuming  perfect  measurements  and  no  noise  this  scheme  works  very 
well.  Figure  3-9  shows  two  M  trajectories  for  a  situation  where  a  step 
change  in  w  occurs.  In  this  case  the  size  of  the  step  is  10%  of  the 

^  ^  I  ^  < 

Larmor  frequency  and  is  applied  along  the  Y  axis  at  the  time  when  M  =  |  M| 

Y.  As  can  be  seen,  the  precession  plane  tips  down  ~  5.7  deg  (0.01  rad) 

with  no  control;  but  is  visually  unaffected  when  a  magnetic  field  is 

applied  to  buck  out  the  rotation.  With  control  applied,  the  vertical 

—  8 

component  of  the  magnetization  reaches  only  a  value  of  -3.11  x  10 
| M ( 0 )  j  ,  and  at  a  sample  rate  of  50  samples  per  Larmor  period,  the  cross 
axis  input  is  nulled  after  only  three  samples. 
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Figure  3-9.  Removal  of  Cross  Axis  Effect  with  3  DOF  Gyro  Scheme 
The  horizontal  precession  plane  is  the  plane  for  zero 
cross  axis  rate,  the^  tilted  precession  plane  is  that 
for  (tiy  «  0.1  yB  Y.  Use  of  3  DOF  gyro  control 
scheme  removes  effects  leaving  the  initial 
horizontal  precession  plane  virtually  unchanged. 

There  are  some  problems  with  this  approach.  One  is  that  the  computation 
involved  will  limit  the  sample  rate;  but  probably  more  crucial  is  the 
susceptibility  to  noise.  Figure  3-10  shows  the  source  of  the  noise 
sensitivity.  Here  the  first  and  third  measurements  are  exact  but  the  second 
is  noisy.  Since  the  algorithm  to  determine  co  fits  a  circular  precession 
plane  to  these  three  points,  we  see  that  the  same  amount  of  noise  causes 
more  problems  for  fast  sample  rates  than  for  slow.  So  the  measurement  noise 
also  places  restrictions  on  the  sample  rate. 
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A 

z 


Figure  3-10.  Effect  of  Measurement  Noise  of  3  DOF  Gyro  Scheme 
The  upper  sketch  shows  the  effdct  of  a  noisy  measurement 
for  a  sample  rate  with  respect  to  precession  rate. 

The  computed  precession  plane  is  orthogonal  to 
the  actual  precession  plane.  The  lower  sketch 
shows  che  effect  of  the  same  amount  of  noise 
when  the  sample  rate  is  slow.  Here  the 
actual  and  computed  precession  planes 
are  nearly  the  same. 


29 


One  final  problem  with  this  scheme  should  also  be  mentioned;  and  that  is 
the  effect  of  magnetometer  drift  on  stability  of  the  applied  field.  Since 
the  SQUID  magnetometers  require  a  small,  but  finite  current  in  the  pickup 
coil,  they  will  also  generate  magnetic  fields.  The  effects  of  the  pickup 
coil  currents,  transverse  to  the  main  field  coil,  will  be  small;  but  the 
pickup  coil  current  in  the  same  axis  will  add  to  the  applied  field 

directly.  If  this  current  were  constant,  its  effect  could  be  modeled  out; 

-9  -5 

however,  a  current  change  of  1x10  amps  gives  a  drift  of  /v  5x10 

rad/sec  for  a  typical  device. 

THREE  ORTHOGONAL  SINGLE  DEGREE-OF- FREEDOM  (SDOF)  GYROS 


We  now  return  to  see  what  we  can  do  with  three  orthogonal  SDOF  gyros, 
each  having  one  or  two  magnetometers  orthogonal  to  each  other  (if  two)  and 
to  the  applied  field  axis.  Since  we  now  have  added  field  coils  in  the  cross 
axes  to  null  cross  axis  rates,  this  is  not  the  same  problem  as  was 
previously  investigated.  These  added  field  coils  take  out  different 
components  of  the  body  rate  for  each  of  the  three  gyros.  So  while  the 
solution  for  to  is  still  the  intersection  of  three  spheres  the  centers  of 
the  three  spheres  are  no  longer  fixed.  Just  as  the  applied  field  moves  the 
center  of  a  given  sphere  along  the  negative  input  axis  for  that  particular 
gyro,  so  too  do  the  feedback  fields  move  the  center  along  the  cross  axes. 

The  question  of  ambiguities  is  not  as  clear  here  as  it  was  for  the  case 
with  no  cross-axis  field  coils.  However,  one  simple  algorithm  appears  to 
work  very  well  for  determining  o>  from  the  three  gyros. 
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The  procedures  is  as  follows: 

1)  Measure  the  precession  frequency  for  each  gyro. 

2)  Compart  is  measured  frequency  with  the  zero  input  Larmor  frequency 
and  take  the  difference  to  be  the  rate  about  the  input  axis. 

3)  Use  these  estimated  rates  to  determine  the  required  cross  axis 
fields. 

4)  Continue  this  procedure. 

Figure  3-11  shows  the  response  to  a  step  in  OJ  .  The  step  size  is  25%  of 

A 

the  Larmor  frequency  and  is  applied  along  the  Z  axis.  The  significance  to 

the  time  scale  on  Figure  3-11  is  that  the  precession  frequency  of  each  gyro 

is  assumed  to  be  determined  in  one  time  unit.  So  let  us  look  at  what  is 

involved  in  determining  the  precession  frequency  from  the  SQUID 

magnetometer (s)  output.  If  one  assumes  that,  during  the  measurement  period, 

the  precession  frequency  is  constant  and  the  gyro  housing  rotation  rates  are 

small  compared  to  the  Larmor  frequency,  then  the  SQUID  output  will  be  a 

sinusoidal  signal,  within  a  certain  band,  centered  at  the  Larmor  frequency. 

If  the  unwanted  low-  and  high-frequency  bands  are  filtered  out  then  the 

SQUID  output  M  (t)  will  have  the  form 
s 

M  (t)  =  A  sin  ojt  +  B  cos  hit .  (3-15) 

s 
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Consider  sampling  this  signal  at  a  fixed  sample  period,  T  sec.  Then  it  is 

shown  that  the  nth  sample  is  given  by 

M  (n)  *  b6  (T)  +  (A  sin  to  T  -  B  cos  to  T)  6  (T-l) 
s 

+  2  cos  coT  M  ( n-1 )  -  M  (n-2).  (3-16) 

s  s 


We  can  write  this  out  in  matrix  form  as 


—  — 

- 

M  (0) 
s 

1 

0 

0 

0 

M  (1) 

S 

0 

1 

M  (  0  ) 
s 

0 

B 

M  (2) 

s 

0 

0 

M  (1) 
s 

M  (0) 
s 

A  sin  toT-B  cos  to  T 

M  (3) 
s 

0 

0 

M  (2) 
s 

M  (1) 
s 

2  cos  to  T 

M  (4) 

s 

0 

0 

V3) 

M  (2) 
s 

• 

* 

• 

• 

• 

, 

-1 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

_ 

The  upper  left  partition  is  due  to  initial  conditions  so  if  we  wait  two 


sample  periods  we  have 


—  — 

M  (2) 

M  (1)  M  (0) 

s 

s  s 

M  (3) 

M  (2)  K  (1) 

s 

s  s 

M  (4) 

M  (3)  M  (2) 

s 

s  s 

•  • 

* 

•  • 

cos 


-1 


(3-18) 


Now  define  2  cos  o)T  *  a^  and  -1  *  and  consider  the  case  where  we 
have  noise  y  such  that 


y  *  pa  +  y  (3-19) 

We  know,  via  measurements,  both  y  and  A}  and  we  can  get  a  best  least 
squares  estimate  of  a  as 

a  T  -1  T 

a  =  (p  p)  p  y  (3-20) 

Since  we  know  a2  =  -1  we  could  consider  the  convergence  of  to  -1  as  a 
simple  test  for  a  good  estimate  of  a^.  Finally,  we  estimate  the 

A 

precession  frequency  from  a^  as 
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A 

A  1  7  ®i 

(a  =  ~  cos  (— )  ( 3-21 ) 

Then  the  estimates  of  to  ,  one  for  each  gyro,  can  be  used  to  determine  the 
fields  necessary  to  cancel  the  cross  axis  rates.  Note  that  if  fields  are 
generated,  there  will  be  changes  in  the  precession  frequencies  of  the  three 
gyros.  Thus  the  next  two  samples  will  again  exhibit  the  transient  terms  in 
Equation  (3-17). 

4.0  A  PRACTICAL  IMPLEMENTATION 

Thu  previous  rh.ipt  **r  Invent  tqatori  tin*  dynnmlon  nf  t  ti«-  ^l!<*  «iy r < i  .uni 
pointed  out  several  problems  due  to  cross  axis  sensitivities.  The  latter 
portion  however  hinted  that  these  effects  could  perhaps  be  removed  with  a 
3-orthogonal  gyro  configuration.  Clearly  this  implementation  would  wotk 
with  the  gyros  mounted  on  a  stable  platform.  This,  however,  would  defeat 
the  purpose  of  developing  a  no-moving  parts  technology  and  would  be 
difficult  to  implement  in  an  all  cryogenic  inertial  measurement  unit.  The 
use  of  additional  applied  fields  to  null  out  cross  axis  rates  is  the 
equivalent  of  a  stable  platform  approach  and  requires  no  moving  parts. 

The  basic  strategy  then  is  to  measure  the  magnetization  in  each  gyro  and 
from  these  measurements  estimate  each  of  the  Larmor  precession  frequencies. 
So  the  first  task  is  to  implement  an  estimator  based  on  a  measurement  from  a 
SQUID  sensing  one  component  of  M.  The  most  general  expression  for  the 
component  of  magnetization  sensed  by  the  SQUID  is 

M  (t)  =  A  sin  cu  t  +  B  cos  cut  +  C  (4-1) 

s 
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Taking  the  z-transform  yields 


Sin  0)T„ 

M  (z)  =  A  - 

l-(2cos  coT)  Z  +Z 


jo  l(cos  CO  T)  Z~ _  _1 _ 

+B  -1  +2  +("  -1 
l-(2cos  Oi  T)Z  +Z  1-Z 


(4-2) 


Then  from  Equation  (4-2)  the  nth  sample  is  given  by 


M  (n)  =  (B+C)  6  ( T )  +  ( Asin  CO  T- (  2C  +  B )  COS  OlT-B)  6  (T-l ) 
s 

+  ( C+Bcos  oj  T-Asin  6j  T)  5  (T-2) 

+  (2  cos  ojT  +  1)  M  ( n-1 )  -  (2  cos&,T+l)M  (n-2)+M  (n-3)  (4-3) 

3  S  S 


In  matrix  form  this  becomes 


M  (0) 

s 

M  (1) 

s 

M  (2) 
s 

•  •  •  • 

M  (3) 

s 

M  (4) 

s 

M  (5) 
s 

M  (6) 

s 


M  (n) 

s 


10  0 


0  10 


0  0  1 


M  (0)  0  0 

s 

M  (1)  M  (0)  0 

s  s 


M  (2)  M  (1)  M  <0) 
s  s  s 

M  (3)  M  (2)  M  (1) 
s  s  s 

M  (4)  M  (3)  M  (2) 
s  s  s 

M  (5)  M  (4)  M  (3) 
s  s  s 


sincoT  -  ( 2C+B) cos  coT  +  B 


+  Bcos  coT  -  Asinoj  T 


M  (n-1)  M  ( n-2 )  M  (n-3) 
s  s  s 


cos  coT  +  1 


( 2  cos  T  +  1 ) 


(4-4) 


Finally  removing  the  transient  terms  gives 


• 


.  *  •  V  *  V 


2  cos  T+l 


Ms(3f 

~M  (2) 
s 

M  ( 1 ) 

3 

«.<or 

Ms<4) 

M  (3) 
s 

M  (2) 
s 

Mg(l) 

M  (5) 
s 

M  (4) 
s 

M  (3) 
s 

M  (2) 
s 

M  (6) 
s 

M  (5) 
s 

M  (4) 

s 

* 

• 

* 

M  (3) 
s 

« 

• 

• 

• 

* 

M  (n) 

s 

M  (n-1 ) 
s 

• 

• 

M  (n-2) 

s 

• 

• 

• 

M  ( n-3 

s 

—  _ 

_ 

-(2  cos  T+l) 


(4-5) 


P 


a 


Note  that  the  addition  of  the  steady  state  term  to  M  (t)  adds  another 

s 

parameter  to  estimate  and  correspondingly  adds  another  column  to  the  P 
matrix.  Again  though  the  best  least  squares  estimate  of  a  is 


A 

a 


.  T  ,-l  T 

(p  p)  p  y 


(4-6) 


Now  a^  should  converge  to  1  and  either  a^  or  a2  can  be  used  to 
estimate  0). 

A  program#  written  in  BASIC,  is  included  in  the  appendix  which  performs 
this  estimation.  The  program  simulates  the  dynamics  of  3-orthogonal  3He 
gyros  via  theoretically  derived  state  transition  matrices  to  describe  the 
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motion  of  the  magnetization  vector.  From  this  the  SQUID  measurements  are 
determined  and  finally  a  is  computed.  For  reasons  not  yet  understood  the 
two  parameter  estimator  Equation  (3-10)  was  sometimes  superior  to  the  three 
parameter  estimator  Equation  (4-6).  Thus  both  estimators  are  exercised  and 
convergence  criteria  for  each  evaluated  to  choose  the  better  estimate  of  to  . 

The  approach  appears  to  offer  some  promise  for  accomplishing  the 
estimation  of  to  with  more  than  adequate  bandwidth  for  most  applications 

A 

requiring  precision  gyros.  The  remainder  of  the  effort  would  be  to  use  to 
to  generate  the  cross  axis  fields  in  each  gyro  to  remove  cross  coupling 
terms.  An  optimum  control  law  would  have  to  be  developed  to  accomplish  to 
this  task.  This  work  was  not  completed  due  to  termination  of  the  project. 

5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  problems  inherent  with  rotation  sensing  via  observation  of  free 

precession  of  a  nuclear  species  seem  largely  solvable  by  mounting  the 

devices  on  a  stable  platform.  Since  the  equivalent  of  a  stable  platform  can 

be  accomplished  by  generating  cross  axis  magnetic  fields  it  is  a  promising 

technology  for  high  precision  rotation  measurement  with  no  moving  parts. 

The  3He  device,  in  particular,  promises  high  resolution  and  used  in 

conjunction  with  other  cryogenic  instruments  (i.e.,  accelerometers,  gravity 

gradiometers,  computers,  and  clocks)  could  provide  an  all  cryogenic  inertia] 

14 

measurement  unit  for  very  demanding  applications. 

There  is  a  fair  amount  of  work  to  be  finished  before  the  promise  becomes 
a  reality  however.  From  this  end  additional  effort  is  required  to  optimize 
the  estimation  scheme  and  then  to  develop  and  optimize  the  control  law. 
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From  the  hardware  point,  the  Stanford  device  should  be  assembled  and  tested 
to  verify  the  expected  sensitivity.  Finally,  a  test  of  the  overall  scheme 
should  be  performed,  substituting  the  Stanford  device,  with  actual 
measurements,  for  one  of  the  simulated  devices. 
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10  DIM  A(3),P(50,3),V(3,3),C(3),T(3,3),Y(50),U(3,3) 

11  DIM  Z(3,3),N(3),N(3,3),G(3,3),C(3,3),N(3,3) 

12  DIM  0(53,3) 

30  FOR  Iai  TO  3 

31  FOR  J»1  TO  3 

32  LET  G(l,J)aO 

33  NEXT  J 

34  NEXT  1 

40  PRINT  "ENTER  THE  NUMBER  CF  ROWS  FOR  THE  P  MATRIX" 

41  INPUT  N 

42  LET  K«3 

50  PRINT  "ENTER  THE  LARMGR  FREQUENCIES  FOR  THE  X , I, AND  Z  GYROS" 

51  INPUT  G ( 1 # 1 ) 

52  INPUT  0(2,2) 

53  INPUT  G(3,3) 

54  PRINT  "ENTER  M(0)  FOR  THE  X  GYRO  X,Y,Z  COMPONENTS" 

55  INPUT  M ( 1 , 1 ) 

56  INPUT  M(2,l) 

57  INPUT  M(3,l) 

58  PRINT  "ENTER  M(0)  FOR  THE  Y  GYRO  X,Y,Z  COMPONENTS" 

59  INPUT  H(l,2) 

60  INPUT  M (2,2) 

61  INPUT  M(3,2) 

62  PRINT  "ENTER  M(0)  FOR  THE  Z  GYRO  X,Y,Z  COMPONENTS" 

63  INPUT  M( 1 ,3) 

64  INPUT  M (2 ,3) 

65  INPUT  M(3,3) 

66  PRINT  "  ENTER  THE  TIME  STEP" 

67  INPUT  T 

70  PRINT  "ENTER  THE  ROTATION  RATE  X,Y,Z  COMPONENTS" 

71  INPUT  W(l) 

72  INPUT  «(2) 

73  INPUT  W(3) 

97  REM  «•*#•  MMimtMlMIUI ••*»!#» II ••#•»•« IMMMMMMIf Mil *#«•«•* 

98  FOR  J*1  TO  (N*3) 

99  REM  ***************  X  GYRC  *************** 

100  LET  C(1,3)BG(1,1)+M(1) 

101  LET  C(2,3)«G(2,1)*W(2) 

102  LET  C(3,3)*G(3,1)+N(3) 

110  GCSUB  500 

120  LET  N(l,l)BZ(l,l)*M(l,l)*Z(l,2)*M(2,l)+Z(i,3)*N(3,l) 

121  LET  N(2,l)sZ(2,l)*M(l,l)*Z(2,2)*H(2,l)*Z(2,3)*M(3,l) 

122  LET  N(3,l)sZ(3,l)*M(l,l)+Z(j,2)*M(2,l)+Z(3,3)*M(3,l) 

125  LET  M ( 1 , 1 )*N ( 1 , 1 ) 

126  LET  N(2,t)«N(2,l) 

127  LET  M(3,1)«N(3,1) 

140  REM  **************  Y  GYRC  *************** 

141  LET  C(1,3)BG(1,2)+W(1) 

142  LET  C(2,3)«G(2,2)+N(2) 

143  LET  C(3,3)bg(3,2)+W(3) 

145  GOSUB  500 

150  LET  N(l,2)aZ(l,l)*M(l,2)*Z(l,2)*M(2,2)*Z(l,3)*M(3,2) 

151  LET  N(2,2)bZ(2,1)*M(1,2)*Z(2,2)*M(2,2)*Z(2,3)*M(3,2) 

152  LET  N(3,2)bZ(3,1)*M(1,2)+Z(3,2)*M(2,2)+Z(3,3)*M(3,2) 

155  LET  M(1,2)«N(1,2) 

156  LET  M ( 2 , 2 )*N ( 2 , 2 ) 

157  LET  M(3,2)«N(3,2) 

180  REM  **************  Z  GYRC  *************** 

181  LET  C(1,3)>G(1,J)+W(1) 

182  LfcT  C(2,3)«G(2,3)+N(2) 

183  LET  f(3,3)«G(3,3)*W(3) 

185  GUSUB  500 

190  LET  N(1,3)BZ(1,1)*M(1,3)+Z(1,2)*W(2,3)+Z(1,3)*M(3,3) 

161  T.ET  Nf9.31«Zf9.1  1*MM  .  T  1  *7  f  9 . 9 1  *M  f  9 . 3  W'/ f  9  .  11*Mf  3.31 


192  LET  K(3,3)sZ»,l)M(l,3)42(3,2)«t(2,3)*2(3,3)»M3,3) 

195  LET  M  (  1 , 3  ) *N (1,3) 

196  Ltl  M(2,3)»N(2,3) 

197  LET  M  (3#  3)*N(3,3) 

200  LET  C(J,1)»M(2#1) 

201  LET  Q(J,2)*M(3,2) 

202  LET  G(J, 3)*M(1,3) 

205  NEXT  J 

210  REM  . . 

250  REM  ******  MOM  READY  TO  ESTIMATE  THE  LARMOR  FREQUENCIES  ****** 

251  REM 

252  REM 

255  REM  ***************  x  GYRO  *************** 

256  REM 

260  FOR  Jsl  TO  N 

261  LET  ¥(J)*<I((J*3),1) 

262  LtT  P(J,t)*0((J*2),li 

263  LET  P(J,2)sU((U+l),l> 

264  LET  P(J,3)*U(J,1) 

270  NEXT  J 

275  GOSUB  967 
280  GOSUB  910 
285  GOSUB  971 
290  GOSUB  890 

295  PRINT  A(1),A(2),A(3),C 

296  GOSUB  650 

297  IF  M*3  GO  TO  350 

298  GOSUB  2000 

350  REM  ***************  y  GYRO  *************** 

360  FOR  Jsl  TO  N 

361  LET  Y ( J)c£ ( l  J*3) ,2) 

362  LET  P(J, l)«0((J+2),2) 

363  LET  P(J,2)*U((J*1),2) 

364  LET  P(J,3)*0(J,2) 

370  NEXT  J 

375  GOSUB  987 
380  GOSUB  910 
385  GOSUB  971 
390  GOSUB  890 

395  PRINT  A(1),A(2),A(3),L 

396  GOSUB  650 

397  If  M*3  GO  TU  450 

398  GUSUB  3000 

450  REM  ***************  i  GYRO  *************** 

460  FOR  J*1  TQ  N 

461  LET  Y(J)s<l((J*3 ),3) 

462  LET  P(J,l)«<J((J  +  2),3) 

463  LET  P(J,2)*0((J+1),3) 

464  LET  P<.J,3)*U(J,3) 

470  NEXT  J 

475  GUSUB  987 
480  GOSUB  910 
*85  GUSUB  971 
490  GOSUB  690 

492  PRINT  M1)«A(?),A(3),D 

493  GOSUB  650 

495  IF  M*3  GO  TO  605 

496  GOSUB  4000 

497  GO  TO  605 

605  INPUT  X 

606  IF  A*0  GO  TU  70 
610  GO  TO  <*999 

640  REM  **************************************** *************9 ********** 

641  REM  *  THIS  SUhROUllHF  TESTS  TO  SFE  If  THE  Tw;j  OH  THE  THREE  * 

642  PEP  *  PARAMETER  FST IHATTR  *ORKS  BETTER  V V  * 

643  PEM  ***********************************************************9**9* 


650  IF  ARSfDXl  .00000E-o5  (It  TO  655 

651  IF  AP3(AC2)+A(1))>.1  GC  TC  655 

652  IF  ABS(A(3)-1)>.1  GU  TC  655 

653  LET  Ms 3 

654  GO  TC  660 

655  LET  Ms2 
660  RETURN 

2000  REM  ***************  X  GIRO  * 

2001  REM  2PE 

2002  FCR  J=1  TO  N 

2003  LET  Y{J)=0((J+2),1) 

2004  LET  F(J,1)=0(CJ+1),1) 

2005  LET  P(u,2)=0(J,l) 

2010  NEXT  J 

2015  GOSUB  987 
2020  GCSUB  1000 
2025  GOSUB  971 
2030  GOSUB  890 
2035  PRINT  A(1),A(2) 

2040  LET  M=3 
2045  RETURN 

3000  REM  ***************  Y  GYRO  < 

3001  REM  2PE 

3002  FOR  J=1  TO  N 

3003  LET  Y (J)=0( (J+2) ,2) 

3004  LET  P(J,l)sQl(04l),2) 

3005  LET  P(J,2)=0(J,2) 

3010  NEXT  J 

3015  GGSUE  987 
3020  GCSUB  1000 
3025  GOSUB  971 
3030  GOSUB  890 
3035  PRINT  A(1),A(2) 

3040  LET  M=3 
3045  RETURN 

4000  REM  ***************  Z  GYRO  < 

4001  REM  2PE 

4002  FCR  Jsl  TO  N 

4003  LET  YCJ)=Oi(J+2),3) 

4004  LET  P(J,l)sO((J+l),3) 

4005  LET  P(J,2)=Q(J,3) 

4010  NEXT  J 

4015  GOSUB  987 
4020  GOSUB  1000 
4025  GCSUB  971 
4030  GCSUB  890 
4035  PRINT  AC1),A(2) 

4040  LET  M=3 
4045  RETURN 
9999  END 
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499  REM  *»*********»»*****************»*>****************»*********** 

500  REM  *  SUBROUTINE  PHI  * 

501  REM  *  THIS  COMPUTES  THE  STATE  TRANSITION  MATRIX  FOR  GYRO  ♦ 

502  REM  ************************************************************ 

503  LET  Q«SOR(C(l,3)*C(L,3)*C(2,3)*C(2,3)*C(3,3)*C(3,3)) 

504  LET  A*COS  (Q*T) 

505  LET  B«5IN(G*T) 

510  LET  Z<1,1)*A4C<1,3)4C(1,3)»C1-A>/(Q*0) 

511  LET  Z(l,2)«C<l,3>*C(2,3)*(i-A)/CO*Q)4C(3,3)*B/Q 

512  LET  7<l,3)«C(l,3)*C(3,3>*(l-A)/(Q*G)«CC2,3>*B/0 

513  LET  Z(2,l)*CU,3)*C<2,3)*(l-n/(Q*Q)-C(3,3)*B/G 

514  LET  Z(2,2)«A4C<2,3)*C(2,3)*(1-A)/(Q*0) 

515  LET  Z(2,3)«C(2,3)*C(3,3)*(l-A>/(Q*Q)4C(lr3)*B/Q 

516  LET  Z(3,1)*C(1,3)*C(3,3)*(1-A>/(G*G)4C(2,3>*B/0 

517  LET  Z<3,2)«C(2,3>*C<3,3)4<1-A)/(Q*0)-C(1,3)*B/G 
516  LET  Z(3,3)«A+C(3,3)*C(3,3)*(i-A)/(G*0) 

520  RE1URN 
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980  REM  ********************************************** 

981  REM  *  SUBROUTINE  PTP  * 

982  REM  *  THIS  SUBROUTINE  COMPUTES  T*  P (TRANSPOSE) *P  * 

984  REM  ********************************************** 
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987  FOR  1*1  TO  3 

988  FOR  J*t  TO  3 

989  LET  T (I , J)*0 

990  NEXT  J 

991  NEXT  I 

992  FOR  1*1  TO  M 

993  FOR  0*1  TO  M 

994  FOR  K*1  TO  N 

995  LET  TCI,J)*T(I,J)+P(K,I)*P(K,J) 

996  NEXT  K 

997  NEXT  0 

998  NEXT  1 

999  RETURN 


01234567690123456789 

01234567890123456789 

01234567890123456789 

01234567890123456789 


44  H3X-11M  V3.2  4* 
44  R8X-11M  V3.2  4* 
44  RSX-lt*  V3.2  4* 
44  RSX-11N  V3.2  4* 


Cil3,lJINVERS  -  NO  PAGE  LIMIT 
FORM  #0  -  NORMAL  HARDWARE  FORMS 
NO  IMPLIED  FORM  FEED 
DP0'tli3,l)INVER8.BA8!l 


till 

11 

n 

333333 

11 

c  c  c  c 

11 

u 

333333 

11 

1111 

mi 

33  33 

ini 

1111 

mi 

33  33 

nn 

11 

u 

33 

n 

11 

n 

33 

n 

11 

ii 

33 

9  9  9  9 

n 

11 

u 

33 

9  9  9  9 

n 

11 

n 

33 

9  9  9  9 

n 

11 

n 

33 

9  9  9  9 

n 

11 

n 

33  33 

9  9 

n 

11 

ii 

33  33 

9  9 

n 

({[( 

min 

mm 

333333 

9  9 

nn 

c  c  c  c 

min 

mm 

333333 

9  9 

nn 

r- 


£ 

min 

NN 

NN 

VV 

vv 

EEEEEEEEEE 

RKRRRRRR 

ssssssss 

mm 

NN 

NN 

vv 

vv 

EEEEEEEEEE 

RRRRRRRR 

SSSSSSSS 

ii 

NN 

NN 

VV 

vv 

EE 

RR 

RR 

$s 

1 

ii 

NN 

NN 

vv 

vv 

EE 

RR 

RR 

55 

5 

ii 

NNNN 

NN 

vv 

vv 

EE 

PR 

RR 

SS 

V 

ii 

NNNN 

NN 

vv 

vv 

EE 

RK 

RR 

ss 

> 

V 

ii 

NN 

NN  NN 

vv 

vv 

EEEEEEEE 

RRRRRRPR 

586888 

V 

ii 

NN 

NN  NN 

vv 

vv 

EFEEEEEE 

RRRRRRRR 

ssssss 

$ 

n 

NN 

NNNN 

vv 

vv 

EE 

RR 

RR 

ss 

B 

ii 

NN 

NNNN 

vv 

vv 

EE 

RR 

RR 

ss 

W 

ii 

NN 

NN 

vv 

vv 

EE 

RR 

RR 

ss 

ii 

NN 

NN 

vv 

vv 

EE 

RR 

RR 

ss 

mm 

NN 

NN 

vv 

EEEEEEEEEE 

RR 

RR 

ssssssss 

mm 

NN 

NN 

vv 

EEEEEEEEEE 

RR 

RR 

ssssssss 

1 

01234567890123456789 

** 

RSX-11M 

V  3 , 2 

01234567890123456789 

** 

RSX-ilM 

V3»2 

M 

01234567890123456789 

4* 

R8X-11M 

V3.2 

** 

•• 

/ 

01234567890123456789 

4* 

RSX^J  M 

V3.2 

4* 

i 


tll3'l)XNVCRS  -  NO  PAGE  UNIT*. 
FORM  #0  -  NORMAL  HARDWARE  FORfi 
NO  IMPLIED  FORM  FLED 
DP0|(U3rlJ  IN  VCRS  •BAR  1 1 


i 


909  REM  M»»MMM**M*MM*MMM*****M*«**M**M*M*I**M*IM* 

910  REM  ♦  SUBROUTINE  INTERS 
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912  REM  *  SUBROUTINE  CALLS  ADJNT  AND  DTRMNT 

913  REM 

914  60SUB  953 

915  GGSUB  930 
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919  LET  V(I,J)«U(I,J)/D 
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931  RE#  *  SUBROUTINE  ADJM  * 

933  RE#  *  THIS  SUBROUTINE  CONFUTES  UMDdCINl  Cl )  * 

933  REM  I********************************************************** 

936  LET  U(1,1)>T(2,2)*T(3,3>-1(3,2>*TC2,3) 

937  LET  U(2,l)«TC3,l)*T{2,3)-l<2,i)*TC3,3) 

939  LET  U<3,1)*T(2,1)*T(3,2)-T<3,1)*T(2,2) 

939  LET  U(l,2)«T<3f2)«T(t,3)-T(l,2>*T<3,3) 

940  LET  U(2«2)«T(l»l)*T(3#3)*T(3fl)*T(ir3) 

941  LET  UC3#2)«T(3,ima#2)-T<l,l>*TC3,2) 

942  LET  U(l,3)*T(l#2)*T(2,3)-T(2,2)*Ta,3) 

943  LET  UC2,3)«T(2,l>*T(l,3)»T<if 1)*T(2,3) 

944  LET  U<3,3)«T(1,1)*TC2,2)-T<2#1)4T<1,2) 

945  RETURN 
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